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Abstract. This article is devoted to define and solve an evolution equation of the 
form dyt = Ay t dt + dX t (yt), where A stands for the Laplace operator on a space of 
the form L p (M. n ), and X is a finite dimensional noisy nonlinearity whose typical form is 



given by Xt(<p) = YliLi X \fid-P)i where each x = [x^ l \ . . . , x^ N ^) is a 7-Holder function 
generating a rough path and each fa is a smooth enough function defined on L p (M. n ). 
The generalization of the usual rough path theory allowing to cope with such kind of 
systems is carefully constructed. 



1. Introduction 

The rough path theory, which was first formulated in the late 90's by Lyons 
and then reworked by various authors [T51 [H], offers a both elegant and efficient 
way of defining integrals driven by rough signal. This pathwise approach enables the 
interpretation and resolution of the standard (rough) differential system 

dy t = a(y t ) dx t , y = a, (1) 

where x is only a Holder process, and also the treatment of less classical (rough) differ- 
ential systems such that the delay equation [28] or the integral Volterra systems 0, HO] • 
In all of those situations, the fractional Brownian motion stands for the most common 
process for which the additional hypotheses required during the construction are actually 
satisfied. 

In the last few years, several authors provided some kind of similar pathwise treatment 
for quasi-linear equations associated to non-bounded operators, that is to say of the 
rather general form 

dy t = Ayt dt + dX t (y t ), t € [0, T] (2) 

where T is a strictly positive constant, A is a non-bounded operator defined on a (dense) 
subspace of some Banach space V and X £ C([0, T] x V; V) is a noise which is irregular 
in time and which evolves in the space of vectorfields acting on the Banach space at 
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stake. Their results apply in particular to some specific partial differential equations 
perturbed by samples of (infinite-dimensional) stochastic processes. 

To our knowledge, two different approaches have been used to tackle the issue of giving 
sense to (|2J): 

• The first one essentially consists in returning to the usual formulation (JTJ) by 
means of tricky transformations of the initial system ([2]). One is then allowed to 
resort to the numerous results established in the standard background of rough 
paths analysis. As far as this general method is concerned, let us quote the 
work of Caruana and Friz [I], Caruana, Friz and Oberhauser [5], as well as the 
promising approach of Teichmann [36] . 

• The second approach is due to the last two authors of the present paper, and is 
based on a formalism which combines (analytical) semigroup theory and rough 
paths methods. This formulation can be seen as a "twisted" version of the 
classical rough path theory. 

Of course, one should also have in mind the huge literature concerning the case of 
evolution equations driven by usual Brownian motion, for which we refer to [5] for 
the infinite dimensional setting and to [7] for the multiparametric framework. In the 
particular case of the stochastic heat equation driven by an infinite dimensional Brownian 
motion, some sharp existence and uniqueness results have been obtained in [HT] in a 
Hilbert space context, and in [3] for Banach valued solutions (closer to the situation we 
shall investigate). In the Young integration context, some recent efforts have also been 
made in order to define solutions to parabolic [27J, [18] or wave type [33] equations. We 
would like to mention also the application of rough path ideas to the solution of dispersive 
equation (both deterministic and stochastic) with low-regularity initial conditions [T6] . 

The present article goes back to the setting we have developed in [19], and proposes 
to fill two gaps left by the latter paper. More specifically, we focus (for sake of clarity) 
on the case of the heat equation in R n with a non-linear fractional perturbation, and 
our aim is to give a reasonable sense and solve the equation 



where A is the Laplacian operator considered on some L p (W n ) space (with p chosen large 
enough and specified later on), namely 



Then the first improvement we propose here consists in considering a rather general noisy 
nonlinearity X evolving in a Holder space C 7 (L p (IR n ); L p (M n )), with 7 < 1/2, instead of 
the polynomial perturbations we had in [Tj5]. A second line of generalization is that we 
also show how to push forward the rough type expansions in the semi-group context, 




(3) 



A : D(A) C L p (R n ) -> L p (R n ). 
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and will be able to get some existence and uniqueness results up to 7 > 1/4, instead of 
7 > 1/3. 

As usual in the stochastic evolution setting, we study equation in its mild form, 
namely: 

yt = S t y + S t _ s dX s (y s ), (4) 
Jo 

where S t : L p (M. n ) — > L p (IR n ) designates the heat semigroup on W l . This being said, and 
before we state an example of the kind of result we have obtained, let us make a few 
remarks on the methodology we have used. 

(a) The main price to pay in order to deal with a general nonlinearity is that we only 
consider a finite dimensional noisy input. Namely, we stick here to a noise generated 
by a 7-Holder path x = (x^\ . . . , x^) and evolving in a finite-dimensonal subspace of 
C(L p (R n ); L p (R n )), which can be written as: 

X t (<p) = 4fM, (5) 
i=i 

with some fixed elements {fi}i=i,...,N of C(L p (W n ); L p (R n )), chosen of the particular form 

for sufficiently smooth functions cr, : K™ xl->l. 

Note that the hypothesis of a finite-dimensional noise is also assumed in [lj or |36j . 
Once again, our aim in [TJ] was to deal with irregular homogeneous noises in space, 
but we were only able to tackle the case of a linear or polynomial dependence on the 
unknown. As far as the form of the nonlinearity is concerned, let us mention that [I] 
deals with a linear case, while the assumptions in [36J can be read in our setting as: one 
is allowed to define an extended function fi(<p) := S-tfi(St<p), which is still a smooth 
enough function of the couple (t, tp). As we shall see, the conditions we ask in the present 
article for /, are much less stringent, and we shall recover partially the results of [36] at 
Section [51 

(b) In order to interpret (Jl]), the reasoning we will resort to is largely inspired by the 
analysis of the standard rough integrals. For this reason, let us recall briefly the main 
features of the theory, as it is presented in [14] : the interpretation of J y s dx s (with 
x a finite-dimensional irregular noise) stems from some kind of dissection of the usual 
Riemann-Lebesgue integral J ydx, when x is a regular driving process. This work of 
dismantling appeals to two recurrent operators acting on spaces of /c-variables functions 
{k > 1): the increment operator 5 and its potential inverse, the sewing map A, the 
existence of which hinges on some specific regularity conditions. If y is a 1-variable 
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function, then 5 is simply denned as (Sy) ts = Vt — Us, while if z ts = f^iut — Du) dx u , then 
(Sz) tus = (Sy)tu(Sx) us . With such notations, one has for instance 



y u dx u =[ / dx u )y s + {y t -y u )dx u = ( / dx u J y s + (5 1 ((8y)(8x))) . 



t / pt \ pt / pt 

dx.. \ v.. + (8~ l ((5v)(5x))) 

ts 

Of course, the latter equality makes only sense once the invertibility of 5 has been 
justified. 

During the process of dissection, it early appears, and this is the basic principles of 
the rough path theory, that in order to give sense to J y s dx s , it suffices to justify the 
existence of the iterated integrals associated to x: x} s = f* dx u , x\ s = f* dx u f u dx v , 
etc., up to an order which is linked to the Holder regularity of x. If x is 7-Holder for 
some 7 > 1/2, then only x 1 is necessary, whereas if 7 G (1/3, 1/2), then x 2 must come 
into the picture. 

Once the integral has been defined, the resolution of the standard system 

(8y) ts = / a(y u ) dx u , y = a, (6) 



where a is a regular function, is quite easy to settle by a fixed-point argument. 

(c) As far as (j4]) is concerned, the presence of the semigroup inside the integral prevents 
us from writing this infinite-dimensional system under the general form ([6]). If y is a 
solution of (jl]) (suppose such a solution exists), its variations are actually governed by 
the equation (let s < t) 



J S t _ u dX u (y u ) 



(8y)ts = yt-y s = S t y - S s y + / [S t - U - S s - U ] dX u (y u ) + 

Jo 

which, owing to the additivity property of the semigroup, reduces to 

(8y)ts = atsVs+ I S t ^ u dX u (y u ), (7) 



where a ts = S t - S — Id. Here occurs the simple idea of replacing 5 with the new coboundary 
operator 5 defined by (Sy) ts = (Sy) ts — dt s y s - Equation ([7j) then takes the more familiar 
form 

(Sy)ts = / S t - u dX u (y u ) , y = ip- (8) 

J s 

In the second section of the article, we will see that the operator S, properly extended 
to act on k- variables functions (k > 1), satisfies properties analogous to 5. In particular, 
the additivity property of S enables to retrieve the cohomology relation 55, which is at 
the core of the most common constructions based on 5. For sake of consistence, we shall 
adapt the notion of regularity of a process to this context: a 1-variable function will be 
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said to be 7-Holder in the sense of 8 if for any s,t, \(8y)t s \ < c \t — s| 7 . It turns out that 
the properties of 8 suggests the possibility of inverting 8 through some operator A, just 
as A inverts 8. This is the topic of Theorem 12.101 which was the starting point of [TU] 
and also the cornerstone of all our present constructions. 

(d) Sections 3 and 4 will then be devoted to the interpretation of the integral appearing 
in (jBJ). To this end, we will proceed as with the standard system ([6]), which means that 
we will suppose at first that X is regular in time and under this hypothesis, we will look 
for a decomposition of the integral in terms of "iterated integrals" depending only on 
X. For some obvious stability reasons, it matters that the dissection mainly appeal to 
the operators 8 and A. 

However, in the course of the reasoning, it will be necessary to control the regularity 
in time of the function u 1— > fi(y u ) as a function of the regularity of y. To do so, 
one can only resort to the tools of standard differential calculus, based on the Taylor 
formula. Unfortunately, those methods can't take our (^-formalism into consideration. 
For instance, it would be futile to search for an equivalent of the rule 

S(fi(y)h(0 = f draKCysiO + ridyUCmSyUO, (9) 
Jo 

which should be expressed in terms of 8 only. This obvious remark obliges us to alternate 
the use of the two operators 8 and 8, but the procedure raises some issues as far as Holder 
regularity is concerned. Indeed, a function which is 7-Holder in the classical sense, that 
is in the sense of 8, is not necessarily 7-Holder in the sense of 8. In such a situation, if 
we refer to the definition of 8 ((8y) ts = (8y) ts — cttsUs), we would like to retrieve \t — s| 7 - 
increments from the operator at s itself. This can be done by letting the fractional 
Sobolev spaces come into play. Namely, wet set B = L p (M n ) and for a G [0,1/2), we 
also write B a , P for the fractional Sobolev space of order a based on B (the definition will 
be elaborated on in Section [2]). One can then resort to the relation (see Section [2]) 

iftp G B a , p , \\a ts ip\\ Bp <c\t- s\ a \\(p\\B a , p - 

Of course, we will have to pay attention to the fact that this time regularity gain occurs 
to the detriment of the spatial regularity. It is also easily conceived that we will require 
B a)P to be an algebra of continuous functions, which explains why we work in some LP 
spaces with p large enough. 

The difficulties evoked by equation (Q are specific to the non-linear case. If the 
vectorfields {fi}i=i,...,N are linear, then we don't need any recourse to the Taylor formula 
and the decomposition can be written thanks to 8 and A only. This particular case 
has been dealt with in [19], as well as the polynomial case, for which we suggested a 
treatment based on trees- indexed integral [T^J [T7] . In our situation, we shall see that the 
landmarks of the construction, that is to say the counterparts of the usual step-2 rough 
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path ( f dx, JJ dx <g> dx), are (morally) some operators acting on B, defined as follows: 
for ip.tp G £>, set 



for i,j = l,...,N, where <p • i/j is the pointwise multiplication operator of tp by ip- 

In a quite natural way, the results established in Section [3] by using development at 
first order only, will be applied to a 7-Holder process x with 7 > 1/2. The considerations 
of Section HI which involve more elaborate developments, will then enable the treatment 
of the case 1/3 < 7 < 1/2. Finally a few words will be said about the case 7 G (1/4, 1/3] 
in Section [TJ and we shall see how the stack of operators allowing the rough path analysis 
grows at order 3. 

It is also crucial to see how our theory applies to concrete situations. To this purpose, 
using an elementary integration by parts argument, we will see in Section [6] that in order 
to define the operators given by (fTOl) and (TTTT) properly, the additional assumptions on x 
reduce to the standard rough-paths hypotheses. In this way, the results of this article can 
be applied to a iV-dimensional fractional Brownian motion x with Hurst index H > 1/4, 
thanks to the previous works of Coutin-Qian [6] or Unterberger [ID]. This also means 
that in the end, the solution to the rough PDE ([3]) is a continuous function of the initial 
condition and 3 , which suggests (as [36] does) that one can also solve the noisy 

heat equation by means of a variant of the classical rough path theory. However, we 
claim that our construction is really well suited for the evolution equation setting, insofar 
that the arguments developed here can be extended naturally to an infinite dimensional 
noise, at the price of some more intricate technical considerations. We plan go back to 
this issue in a further publication. 

With all these consideration in mind, we can now give an example of the kind of result 
which shall be obtained in the sequel of the paper (given here in the first non trivial 
rough case for X, that is a Holder continuity exponent 1/3 < 7 < 1/2): 

Theorem 1.1. Let X be a noisy nonlinearity of the form |3J) ; where: 

(i) The noisy part x is a N -dimensional Holder- continuous signal in C 7 ([0, T}; IR ) for 
a given 7 > 1/3. We also assume that x allows to define a Levy area x 2 in the sense 
given by Hypothesis^ 

(ii) Each nonlinearity fi can be written as [fi(<f)](C) — where the function 
Oi : R n x IR — > K is such that ai(-,rj) = outside of a ball B k n(0, M), independently of 
rj G IR. We also ask 77 > rj) to be a Cf (IR) function for all rj G IR^. 




(10) 




(11) 
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Then equation admits a unique solution y on an interval [0, T], for a strictly 
positive time T which depends on x and x 2 . Furthemore, the solution y is a continuous 
function of (yo, x,x 2 ). 

Notice that this theorem is directly applicable to the fractional Brownian setting for 
H > 1/3. The case of a Holder coefficient l/4<7<l/3is also discussed at the end of 
the article. 

Here is how our paper is structured: Section [2] is devoted to recall somme basic facts 
about algebraic integration with respect to a semi-group of operators, taken mainly 
from [19]. Then we deal with the easy case of Young integration at Section [31 This 
allows to solve equations for a noisy input with any Holder continuity exponent 7 > 1/2 
(recall that we had to consider 7 > 5/6 in [19]), and it should also be mentioned that 
we obtain a global solution for the RPDE (jlj) in this case. The first rough case, that 
is a Holder continuity exponent 7 G (1/3,1/2], is handled at Section HI Observe that 
the abstract results obtained there are expressed in terms of the operators X x ,X xa and 
X xx defined at equation ffTUj) and (TTTT) . It is also important to notice that only local 
solutions are obtained in the general case, due to the fact that our nonlinearity cannot 
be considered as a bounded function on the Sobolev spaces B Q ^ P . We will show however 
at Section [5] that considering a smoothed version of the nonlinearity, a global solution 
to equation (jlj) can be constructed. Section [6] is then devoted to the translation of these 
results in terms of x 1 and x 2 by a simple integration by parts argument, and thus to the 
application of the abstract results to concrete situations. Finally, we discuss in Section [7] 
the rougher case of a Holder continuity exponent of the noise x satisfying 1/4 < 7 < 1/3. 



2. Algebraic integration associated to the heat semigroup 

This first section aims at introducing the framework of our study, as well as the 
different tools evoked in the introduction. The main point here is the definition and the 
basic properties of the infinite-dimensional coboundary operator <5 already alluded to in 
the introduction, together with its inverse A. At first, we will recall some elementary 
properties of the heat semigroup, which will actually be used in the construction of A 
(Theorem [2/10]). 

2.1. Framework. We will focus on the case of the heat equation on L p (M. n ), for some 
integrer p that will be precised during the study. We denote by A = A p the Laplacian 
operator, considered on the (classical) Sobolev space W 2 ' p {R n ), and by S t the associated 
heat semigroup, which is also defined by the convolution 

S t <p = g t * V , with g t (0 = _L^H?l 2 /2i. (12) 
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As explained at point (d) of the introduction, the interplay between the linear and 
the non-lienar part of the equation invites us to let the fractional Sobolev spaces come 
into the picture: 

Notation 2.1. For any a > 0, for any p G N* ; we will denote by B 0)P the space 
(Id—A)~ a (L p (M. n )), endowed with the norm 

IMIb q ,p = IMUp(R") + II( _ ^) q v ? IIlp(]r™)- 

Set also B p = B 0:P = L p (R n ) for any p G N* U {oo}. 

The space B a , p is also refered to as the Bessel potential of order (2a, p). Adams ([!]) 
or Stein ([31]) gave a thorough description of those fractional Sobolev spaces. Let us 
indicate here the two properties that we will resort to in the applications: 

• Sobolev inclusions: If < \i < 2a — then B a ^ p is continuously included in the 
space C 0, ^{R n ) of the bounded, /i- Holder functions. 

• Algebra: If 2ap > n, then B a:P is a Banach algebra, or in other words Hv^IIbc,,,, < 

1 1 1 1 $a ,p 1 1 1 1 ,p ' 

The general theory of fractional powers of operators then provides us with sharp 
estimates for the semigroup S t (see for instance [32] or [TT]): 

Proposition 2.2. Fix a time T > 0. St satisfies the following properties: 

• Contraction: For all t > 0, a > 0, St is a contraction operator on B a>p . 

• Regularization: For all t G (0, T\, a > 0, S t sends B p on B a ^ p and 

II^||b QiP < c a , T r a y\\ Bp . (13) 

• Holder regularity. For all t G (0, T\, cp G B a ^ p , 

\\S*P - Vh p < c a) Tt a y\\ Bci}P - (14) 

\\&sM\B P <c a , T t- 1+a y\\ Bci , p . (is) 

At some point of our study, the interpretation of the integral J S tu dx l u fi(y u ) will 
require a Taylor expansion of the (regular) function /j. As a result, pointwise multipli- 
cations of elements of B p are to appear, giving birth to elements of B p /k, k G {1, . . . ,p}. 
In order to go back to the base space B p , we shall resort to the following additional 
properties of S t , which accounts for our use of the spaces B p (p > 2) instead of the 
classical Hilbert space B2: 

Proposition 2.3. For all t > 0, k G {1, . . . ,p}, (f G B p /k, one has 

n(k-l) 

\\St<p\\t3 p < C k ,nt 2 p \\<p\\B p/h , (16) 

\\AS t ip\\ Bp < c k , n t *p \\ip\\B p/k - (17) 
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Proof. Those are direct consequences of the Riesz-Thorin theorem. Indeed, for any 

f e Bp/ k , 

n(fc-l) 

\\St<p\\B P < \\gt *v\\b p < \\gth pnp _ k+1) hh p/k < c k>n t * |M| Sp/fc . 

The second inequality can be proved in the same way, since ASt^p = (^f)'-P — dtgt * 
<p. □ 

Let us finally point out the following result of Strichartz [35], which will be at the 
core of our fixed-point argument through Proposition 13.41 (see also [20] for more general 
results): 

Proposition 2.4. For all a E (0, 1/2), for all p > 1, set 

2 \ V2 



= ( fr- l - Aa f / |/(£ + nj) - /(0| *7 

WO L^|r;|<l 



T/ien / G i3 a>p i/ and only if f £ B p and T a f G S p , and 

ll/llft w ~ll/lk + l|r«/lk. 

2.2. The twisted coboundary 5. Notice that we shall work on n th dimensional sim- 
plexes of [0, T], which will be denoted by 

S% = {(s 1} . . . , s n ) G [0, T] n ; Sl <s 2 <---s n }. 

We will also set C n = C n (S^, V) for the continuous n-variables functions from S% to V, 
for a given vector space V. Observe that we work on those simplexes just because the 
operator S t - U is defined for t > u (i.e. on <Sf.) only. 

Let us recall now two basic notations of usual algebraic integration, as explained in 
[T4] and also recalled in [TH]: we define first an coboundary operator, denoted by 5, 
which acts on the set C n = C n (S^, V) of the continuous n-variables functions according 
to the formula: 

n+l 

5:C n -> C n+l , (5g) tl .., n+1 = ^( -1 )'//;,., ..,., (18) 

i=l 

where the notation ti means that this particular index is omitted. In this definition, V 
stands for any vector space. Next, a convention for products of elements of C n is needed, 
and it is recalled in the following notation: 

Notation 2.5. If g G C n (C(V,W)) and h G C m (W), then the product gh G C m+n -i(W) 
is defined by the formula 



[gh)t 1 ...t m+n - 1 — gti...t„ht 



n • n-\-m — 1 * 
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In point (b) of the introduction, we (briefly) explains why the standard increment 8 
was not really well-suited to the study of (j4j). We will rather use a twisted version of 8, 
denoted by 8, and defined by: 

Definition 2.6. For any n G W,y G C n (B a ,p), for allt\< . . . < t n+ i, 

(Sy)tn+i...ti = (<fyk+i...t! - a tn+ltn y tn ... tl , with a ts = S t _ s -Id sis <t. (19) 

The operator a : (t, s) i— > a ts is only defined on the simplex {t > s}. As a consequence, 
we will have to pay attention to the decreasing order of the time variables throughout 
our calculations below. Note that we will often resort to the notation St 8 for St- S , so as 
to get a consistent notational convention for the indexes. 

The rest of this subsection is devoted to the inventory of some of those results. The 
associated proofs can be found in [T9] . 

Let us start with the fundamental property: 

Proposition 2.7. The operator 8 satisfies the cohomolgical relation 88 — 0. Besides, 
Ker8 ]Cn+1 (B a>p ) = Im8 lCn (B a , P )- 

Now, let us turn to a more trivial result, which will be exploited in the sequel. Re- 
member that we use the notational convention 12.51 for time variables. 

Proposition 2.8. If L e C n -x{V) and M e C 2 {C{V)), then 

8{ML) = (8M)L - M(8L). (20) 

The following result is the equivalent of Chasles relation in the 8 setting. It is an 
obvious consequence of the multiplicative property of S. 

Proposition 2.9. Let x a differentiate process. If yts = J St u dx u f u , then (8y)tus — 
for all s < u < t. 

From an analytical point of view, the notion of Holder-regularity of a process should 
be adapted to this context, and thus, we define, for any a G [0, 1/2), p G N*, K G (0, 1), 

C^{B a , p ) := {y G dCBo,) : sup ^hl^A < (2 l) 

s<t \t — S\ 

Let us take profit of this subsection to introduce the Holder spaces commonly used in 
the fc-increment theory. They are the subspaces of Ci(V), C 2 (\^) and C 3 (V) respectively 
induced by the norms (V stands for any normed vector space): 
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N[y;C»(V)] := inf |£ W\ : y = XYj- 

Now, let us state the main result of this section, which allows to invert the twisted 
coboundary operator 5. 

Theorem 2.10. Fix a time T > 0, a parameter k > and let fx > 1. For any 

h G C3 ([0, T]; B K)P ) D ifer5|c 3 (B ), i/iere exzsfo a unique element 

Ah e n aeM C^ a ([0,T];B K+atP ) 

such that 5 (Ah) = h. Moreover, Ah satisfies the following contraction property: for all 
a G [0, fx), 

M[Ah;Cr a ([0,T];B K+ajP )} < c ajlMiT N[h;C%([0,T};B K , p )}. (22) 

The link between the operator A and a more classical formulation of the rough inte- 
gration theory by means of Riemann sums, is given by the following result: 

Proposition 2.11. Let g G C2(B titP ) such that 5g G Ct^(B KtP ) with fi > 1. Then the 
increment 5f = (Id—A5)g G C2(B KtP ) satisfies 

(5f) ts = lim S tt k+1 9t k+1 t k mB K , p , 

\LLts ~ *U 

(tk)eu ts 

for all s <t. 

3. Young case 

The appellation 'Young case' is used in this section order to indicate that only expan- 
sions at first order will be involved in this section. Although this kind of considerations 
has already been explored in [18] under more general hypotheses concerning the spatial 
regularity of the noise, we think that it is worthwhile to illustrate in a simple setting 
the adaptation of the dissection method to the convolutional context. We will see in 
Section [6] that the general result of Theorem 13.61 can be applied to a noise generated by 
a (finite-dimensional) 7-Holder process x, with 7 > 1/2. This is an improvement with 
respect to [19], where the unnatural condition 7 > 5/6 had to be assumed. 

Throughout this section, we fix a parameter 7 G (1/2, 1), which (morally) represents 
the Holder regularity of x. 

3.1. Interpretation of the integral. The aim here is to give an interpretation of the 
twisted Young integral f S tu dx u z u in terms of 5 and A, and to do so, we shall follow 
the same reasoning as in [HJ [19]: we assume first that x and z are smooth processes, 
and obtain a dissection of the integral S tu dx u z u in terms of 5 and A in this particular 
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case. This allows then to extend the notion of twisted integral to Holder continuous 
signals with Holder continuity coefficient greater than 1/2. 

Thus, assume, at first, that x is real valued and regular (for instance lipschitz, or even 
differentiable) in time, as well as the integrand z, and look at the decomposition 

S tu dx u z u = (J^ S tu dx^j z s + S tu dx u (5z) us . (23) 
Now, if we set r ts = J s ' S tv dx v (Sz) vs , one has 

ft f't ru 

(fi'r)tus = I Sf v dx v (5z) vs / Stv dx v {5z) vu Stu I S uv dx v [Sz) vs , 

J s J u J s 

which, using the fact that S tu S uv = S tv , reduces to 

(Sr) tus =(J S tv dx v ^j (5z) us . (24) 

This first elementary step lets already emerge the object which plays the role of the a 
priori first order increment associated to the heat equation, namely 

~t 

Xtl = / Stv dx... 



L ts — / °tv 
J s 

We are then in position to invert 5 in f!24p thanks to Theorem 12. 101 Indeed, one easily 
deduces, owing to the regularity of x and z, 

X x (5z) e C%(B a>p ) for some a e [0, 1/2). 

As a result, we get 

f S tu dx u z u = X?/zi + A ts (X x >* 5z l ) . (25) 

J s 

As in the standard case algebraic integration setting in the Young setting, we now 
wonder if the right-hand-side of (1231) remains well-defined in a less regular context: 

• From an analytical point of view. The regularity assumption of Theorem 12.101 
imposes the condition: for all i £ {1, ... , N}, 

r ,! fe ! £ C^{B^ P ) with a £ [0, 1/2) and /x > 1. 

Therefore, we shall be led to suppose that z l is ^-Holder (in the classical sense), 
with values in a space B a i )P to be precised, or in other words z % £ Cf(z5 Q / jP ), and 
we will also assume that X x ' % £ C\ (C(B a i tP , B a>p )), with k + 7 > 1. In fact, we 
will see that changing space is not necessary when we apply X x > 1 , so that it will 
be sufficient to consider the case a = a'. 
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• From an algebraic point of view. We know that A is denned on the spaces 
C^(B a)P ) fl Ker 5. This constrains us to assume that 5(X x,t 5z l ) = 0, which, by 
()20p . is satisfied once we admit that 5X X ' 1 = 0. 

Let us record those two conditions under the abstract hypothesis: 

Hypothesis 1. From x, one can build processes X x,z (i £ {1, . . . , N} ) of two variables 
such that, for all i: 

• For any a £ [0, 1/2) such that lap > 1, X x,i £ C^(C(B a)P , B a , P )) 

• The algebraic relation 5X X > 1 = is satisfied. 

Remark 3.1. Actually, the additional condition 2ap > 1 could have been skipped in 
the latter hypothesis. We have notified it so that Hypothesis [TJ meets the more general 
Hypothesis [3] of Section HI 

We are then allowed to use the expression (1251) for irregular integrands: 

Proposition 3.2. Under the assumption (QJ), we define, for all processes z such that 
z % £ Ci(B KjP ) fl Cf (Bp), i = 1, . . . , N , with k < 7 and k + 7 > 1, the integral 

J ts (dx z) = X x t fz i s + A is (X x * 82?) . (26) 

In that case: 

• J{dxz) is well-defined and there exists an element z £ Cj(B KyP ) such that 5z is 
equal to J(dx z). 

• It holds that 

N[z-C1{B^ P )\ < c x {M[z-Cl(B KtP )}+ Miz-^B,)}} , (27) 

with 

c x < c{Af[X x ;C2(£(B p ,B p ))] +N[X x -Cl(C{B K ^B^ p )}} (28) 

• The integral can be written as 

JUSxz) = hm o S ttk+1 X^ ltk zl k , (29) 

(**)6A 

where the limit is taken over partitions Ar s ^i of the interval [s,t], as their mesh 
tends to 0. Hence it coincides with the Young type integral S tu dx u z u . 

Proof. The fact that J ts (dx z) is well defined is a direct consequence of Hypothesis [U 
and the Chasles relation 5J(dxz) 1 which accounts for the existence of z, can be shown 
by straightforward computations using (12"0|) . 
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For the second point, notice that, thanks to Hypothesis [TJ one has 

N[J{dxz)-Cl{B K)P )\ 

< M[X^ C2(C(B K>P , B KtP ))\ A^; C^B^)} + ^[A(X^ Sz*); C 2 7 (i5 K , p )], 

since X x ' 1 5z l G C^ +K (B P ). By the contraction property (122]) of A, we then deduce 

which completes the proof of ( 1271) . According to Proposition 12 .111 (1291) is a consequence 
of the reformultion 

J(dxz) = (Id - A8)(X x ' i z i ). (30) 

□ 

Remark 3.3. The careful readers may wonder if the starting decomposition (|23|) is really 
the most relevant choice as far as the stability of the 5 structure is concerned. Indeed, 
at first glance, it seems more appropriate to look for an expression written in terms of 5 
only, and thus, for a one-dimensional noise x, we would rather rely on the decomposition 

St u dx u z u = StsSxts z s + / St u dx u (5z) us . (31) 

J s 

The order-one operator then becomes Xf s = S ts 5x ts = /* S tu dx u S us , which coincides 
with the first order operator built in [T9] . 

However, one must keep in mind the particular form of the integrand we are about to 
consider in the system (j2j), namely z = f(y), for some non- linear function /. In order 
to settle a fixed-point argument, we will have to control the regularity of this integrand 
according to the regularity of y and to do so, we can only resort to the standard tools 
of differential calculus, which are not consistent with the 5 formalism. In other words, 
if one wants to estimate the norm of Sf(y), one is forced to estimate the norm of the 
classical increment Sf(y) first. In this context, the decompositions (!23l and (I3TT) give 
rise to similar treatments. 

3.2. Resolution of the associated differential system. Using the formalism we have 
just introduced, we are going to show the following result of existence and uniqueness 
of a global solution. To begin with, let us state the assumption on the regularity of the 
functions <7j appearing in the definition of the vectorfields /j, % = 1, . . . , N. 

Hypothesis 2. Let f : B p — > B v be a vector field defined by f(<p){£) = f or 
some function a : lR n xR^I. We say that f e Xk for k > 1 if o is of uniformly 
compact support in the first variable, in the following sense: o~i : W l x R — > R is such 
that o~i(-,rj) = outside of a ball B r n(0, M) , independently of rj G R. 
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In order to be element of X k , we also ask to a vector field f to satisfy the following 
inequality: 

sup max |V™<7j(£, 77)1 + max | V € V™<7j(£, 77) | < +00. 

£eK n ,?7e]R n=0 >"-' fc n=0,...,k-\ 

A direct application of Proposition 12.41 easily leads to: 

Proposition 3.4. If / £ X\, then for any ip £ B a>p , f(<p) £ B a , P and 

AO/%); B a , p ] < c f {1 + M[ip\ B K J} . 

The following notation will also be used in the sequel of the paper. 

Notation 3.5. Let A, B be two positive quantities, and a a parameter lying in a certain 
vector space V . We say that A < a B if there exists a positive constant c a depending 
on a such that A < c a B. When we don't want to specify the dependence on a, we just 
write A < B. Notice also that the value of the constants c or c a in our computations 
can change from line to line, throughout the paper. 

We are now ready to prove the main theorem of this section: 

Theorem 3.6. Assume HypothesisUlwith 7 > 1/2, and assume also that f = (/1, . . . , /at) 
with fi £ X 2 for i — 1, . . . , N. For any n < 7 such that 7 + n > 1 and 2np > n, consider 
the space C°' K ([0, T], B KjP ) = C°([0, T], B KjP ) fl ([0, T], B K ^ P ), provided with the norm 

■A/].; CT' K ([0, T], B KiP )} = Ar[- (*([0, T], B K , P )} + Af[.; Cf([0, T\, B KiP )]. 

Then the infinite-dimensional system 

(k)ts = Jts(dxf(y)) , y = ipeB KjP , (32) 

interpreted with Proposition \3.B, admits a unique global solution in C^QO, T], B K}P ). 
Besides, the ltd application X x ' 1 ) \—> y, where y is the unique solution of [3^1 , is 
Lipschitz. 

Remark 3.7. In the last statement, we consider the operators X x ' 1 as elements of the 
incremental space C1(C{B P , B p )) flC^ (£(B K , p , B K)P )). The regularity of the Ito application 
with respect to X x ' 1 is then relative to the norm 

M\.;CC^} = Af[.; C](C(B P , B P ))\ + Af[.;Cj(jC(B K , p ,B K , p ))}. 

Proof. It is a classical fixed-point argument. We will only prove the existence and 
uniqueness of a local solution. The reasoning which enables to extend the local solution 
into a solution on the whole interval [0,T] is standard; some details about the general 
procedure can be found in [T3| (in a slightly different context). 
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We consider an interval / = [0, T*] with T* a time that may change during the proof, 
and the application T : B KtP ) -> 0^(1, B K)P ) defined by T(y) = and (5T(y)) ts = 

Jts(dxf(y)). 

Invariance of a ball. Let y e C®'^(I, and z = T(y). By ( 1271) . we know that 

Af[z;C?(I,B KjP )\ < c x \I\'<- K {M\fi(y);C?(I,B p )] + M[f i (y);C 1 (B K , p )}} . (33) 

Recalling our convention in Notation 13.51 the assumption /j e Ai is enough to guaran- 
tee that the following bounds holds for ff J\f[fi(<p) — fi(ip);B p ] </ A/"[<£ — ^; B p ] and 
N[fi(tp) ; £> KjP ] </ 1 + A/"[</?; for arbitrary test functions <f,ip. So we have 

WM;Ci(I,Bp)] < f M[y;C^I,B p )] 

< f M[y; C?( J, B^,)] + Af[y; C*(I, B K , P )\ 
< f N[y;Cl'\I^«,p)l 

where, to get the second inequality, we have used the property (IT4"|) of the semigroup. 
We get also U[fi(y);C^(B KiP )} < f 1 + A%; Cl(B KjP )], which, going back to ([33]), leads to 

AA[z;Cf(J,S KiP )] S,/ l/r^ll+^C^J,^)]}. 

Besides, z s = (5z) s0 + S^, hence, since S s is a contraction operator on B K>P , 
N[z-C°{I,B K>P )\ < \I\ K M[z;C?(I,B KtP )] + U\\ Bk , p . 

Finally, 

U[z;C°{ K (I,^, P )} < \\4>h K , P + c x \ir K {l+^C^il^B^)}} . 
Then we choose / = [0, Ti] such that c r T 1 7-K < | to get the invariance by T of the balls 

b*^ = {ye cT' K ([o,r ],^ )P ) : = V>, A%;C^([0,T ],i3 KiP )] < i?}, 
for any T < T 1; with (for instance) R — 1 + 2||'0|| J g (5 . 

Contraction property. Let y e B K , P ) and z = r(y), 5 = T(y). By ( 1271) . 

Cx \ir K {Af[m - m-A{B K , p )]+N[m - mwiM) . (34) 



NON-LINEAR ROUGH HEAT EQUATIONS 17 

In order to estimate the Holder norm Af[fi(y) — fi(y); C*(B P )], we rely on the decompo- 
sition 

VtiO) - <n(Z, UO) ~ <£, i/.(0) + y a (0) 

= *(j/-2/)*.(0 / dr^j/^O + r^CO) 

+ (<^(0 / rfr fo'fc ?/ a (0 + r(5j/)t.(0) - y-(0 + • 

Therefore, 
W,(y) -/,($; crfo)] 

< C/ { A% - y; C°{ K (B KiP )] + A/"[y; C^(^, p )] JV[y - y; C?(S„o)] 



Remember that we have assumed that 2k£> > n, so that, by the Sobolev continuous 
inclusion B K)P C B^ , A% - y; C^B^)] < Af[y - y; C^(B K>P )} and as a result 



WM ~ fi(y);^(B p )} < cAf[ y - y;C°^(B K , p )} [1 + M[y;C° 1 ' K (B K>p )} 
The same kind of argument easily leads to 

M[h{y) - m-A{B K , p )] 

< c A% - y; C°{ K (B K , P )] { 1 + Af[y; C°{\B K , P )} + Af[y; C°{ K (B K , P )] 



The last two estimations, together with provide a control of Af[z — Z]C1{B KtP )\ in 
terms of y, y. Moreover, as z$ = zq = if), 

M[z - z-C^B^)] < \I\ K N[z-~z;C«(B K , p )]. 

Now, if y, y both belong to one of the invariant balls B^ o ^, with To < Ti, the above 
results give 

M[z - z; C°'*([0, T ], B K , P )} < c x Tl K {1 + 2R}M[y - y; ^([0, T ], 

It only remains to pick T < T\ such that c x T^ K {1 + 2i?} < |, and we get the contrac- 
tion property of the application Y : B^ Q ^ — > B^^- This statement obviously completes 
the proof of the existence and uniqueness of a solution to (152"]) defined on [0, To]. 

□ 

4. Rough case 

The aim of this section is to go one step further in the rough path procedure: We 
would like to conceive more sophisticated developments of the integral so as to cope with 
a 7-Holder driving process, with 7 G (1/3, 1/2). 
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4.1. Heuristic considerations. The strategy to give a (reasonable) sense to the inte- 
gral J* S tu dx l u fi(y u ) w iH be largely inspired by the reasoning followed for the standard 
integral f*y u dx u , explained in [TH [19]. Thus, let us suppose at first that the process x 
is differentiable (in time), as a function with values in a Banach space. The procedure 
to reach a suitable decomposition of the integral divides into two steps: 

• Identify the space Q of controlled processes which will accomodate the solution 
of the system. 

• Decompose f* S tu dx l u fi(y u ) as an element of Q when y belongs itself to Q, until 
we get an expression likely to remain meaningful if x is less regular. 

This heuristic reasoning essentially aims at identifying the algebraic structures which 
will come into play. The details concerning the analytical conditions will be checked a 
posteriori. The noisy nonlinearity is given by equation (jSJ), namely 

N 

i=i 

and we shall see that a, has to be considered as an element of X 2 , as defined in Hypoth- 
esis EJ 

Step 1: Identification of the controlled processes. The first elementary decomposition 
still consists in: 

S tu dxif t .(y u ) =( [ S tu dxi) fi(y s ) + / S^Miy))^ (35) 



It is then natural to think that the potential solution of the system is to belong to a 
space structured by the relation 

(6y) u =( I S tu dx{] yf + yt 



with admitting a Holder regularity twice higher than y. For the solution itself, we 
would have yf l = fi(y s ), y\ s = J* S tu dx l u S(fi(y)) us hence the potential algebraic struc- 
ture of the controlled processes 

Q = {y : 5y = Xgtf* + y\ 8 }, with = f S tu dx\. 

J s 

Remember that the latter operator satisfies the algebraic relation 

5X X ' { = 0. (36) 
Besides, it will turn out useful in the sequel to write X x ' % as 

X x,i = X ax,i + 6x i ^ with X ax,i = J ^ d ^ ^ 
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Morally, X ax,% admits a higher Holder regularity than x owing to the property (fl4j) 
of the semigroup. We will go back over the usefulness of this trivial decomposition in 
Remark 14.31 In the following we will omit sometimes the vector indexes i, j, . . . whenever 
the contractions are obvious. 

Step 2: Decomposition of f S tu dx u fi(y u ) when y E Q. Going back to expression (1351) . 
we see that it is more exactly the integral S tu dx u S(fi(y)) us that remains to be 
dissected when y e Q, that is to say when the (^-increment of y can be written as 
($y)ts — Xt/y* ,z + y\ s . To this purpose, let us introduce a new notation which will 
appear in many of our future computations: 

Notation 4.1. For any f G X 2 as defined in Hypothesis^ we set 

[fMKO = V 2 a(^(0), 

where V2 stands for the derivative with respect to the second variable. The function f 
is understood as a mapping from B p to B p for any p > 1. 

Using this notational convention, notice that 

S(fi(y)) ts = (Sy) u ■ fiiVs) + [ dr [f>(y s + r(Sy) ts ) - ■ (Sy) u 

Jo 

= (auy.) ■ fiiVs) + (Sy)ts ■ fi{y s ) + fi(y)ts (38) 
= (a ts y s ) ■ fXv.) + {X?h x s ' j ) ■ fM + h{vt + mi's 
= (auy.) ■ ti(y s ) + (Sx% ■ y^ ■ f[(y s ) + f t (yf£ + f t (yf t s + mi's, 
where we have successively introduced the notations 

Mvt 1 = f dr [f>{y s + r{8y) ts ) - f t (y a )] ■ (8y) ts , /,(„){? = y\ s ■ fl(y s ), (39) 
Jo 

fi(y)l? = (Xu' j ys' j )-my s ). (40) 

Observe that, in the course of those computations, we have used some additional 
conventions that we make explicit for further use: 

Notation 4.2. Let ip,ip be two elements of B p . Then ip-i/; is the element of B p / 2 defined 
by the pointwise multiplication [tp ■ ^>](£) = </ 7 (0^'(6)- If we assume furthermore that 
M,N are two elements of C(B P ;B P ), then the bilinear form B(M <S> N) is defined as: 

B(M ®N):B p xB p ^ Bp /2 , (<p, rj)) i-> [B(M ® N)](<p, V>) = M(<p) ■ N{i(>). 

With this convention in mind, the algebraic decomposition (1381) of fi(y) can now be 
read as: 

6(fi(y)) u = B(a ts ® M)(y, fl(y)) s + (8x% ■ y^ ■ f[(y s ) + h(y)\ s . (41) 
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If we analyze the regularity of the terms of this expression, it seems reasonable to 
consider the first two terms as elements of order one and fi(y)^ as an element of order 
two. Let us make two comments about this intuition: 

(a) To assert that B(a ts <g> ld)(y, fl(y)) s admits a strictly positive Holder regularity, 
otherwise stated to retrieve increments \t — s\ a from the operator at s , we must 
use the property (114j) of the semigroup. It means in particular that a change of 
space will occur: if y s G £> QjP , then B(a ts ® ld)(y, fl(y)) s will be estimated as an 
element of B p . This remark also holds for fi(y)lf = (X^^y^) ■ f-(y s )- 

(b) The term fi(y)^' 1 is considered as a second order element insofar as it is easily 
(pointwise) estimated by (a constant times) \(5y)t s \ 2 - However, as far as the 
spatial regularity is concerned, this supposes that /i(y)"' 1 has to be seen as an 
element of B p / 2 , if y G B p . To go back to the base space B p , we shall use 
the regularization properties (fT6l) of the semigroup, through the operator X x 
(Hypothesis (1531). 

Now, inject the decomposition (|4"T|) into fl35|) to obtain 

f Stu dxi f t (y u ) = X-'Myl + X?\y t f>(y)) s + X^\y^ • //(„)), 

+ t StudxiMyt,, (42) 

J s 

where we have introduced the following operators of order two (which act on some spaces 
which will be detailed later on): 

X™> % = f S tu dxiB(a us ® Id) and X™> ij = f S tu dx\ (5rf) us . (43) 

J s J s 

A little more specifically, those operators act on couples (if, ip) in some Sobolev type 
spaces, and 

X? s a, \(p,ip) = / S tu dx\ [a us (ip) ■ ip] and X^ (<p) = / S tu dx\ (5x J ) us [ip\ . 

J s J s 

Then, since we have assumed that fi(y)^ admitted a "double" regularity, we can see 
the residual term r ts = J s S tu dx l u fi(y third order element, whose regularity is 

expected to be greater than 1 as soon as the Holder regularity of x is greater than 1/3. 
Thus, we are in the same position as in (123]) . and just as in the latter situation, r will 
be interpreted thanks to A. 

In order to compute 8r, rewrite r using (j42p : 

rts = f s S tu dx\ U{Vu) - XgfM - X^\y, f'{y)) s - X^\y^ • f[{y)) s . 
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Therefore, with the help of the algebraic formula (1201) . we get 

(5r) tus = XtMiy))^ - (5X* a '% us (y, f'(y)) s + X^5(y, f'(y)) us 

- (5X™'% us (y*'i ■ ti(y)) s + X^ ij 5(y*j ■ f>(y)) us . 

Going back to the very definition of X xa ' 1 and X xx ' % \ it is quite easy to show that the 
following relations are satisfied whenever x is a smooth function: 

(6X* a '% us = X™>\a us ® Id) + X^(a us ® Id), (44) 

{8X™>% US = X*£{8x?) u8 . (45) 
By combining these two relations together with ([41]) . we deduce 

(Sr)tus = X^iml) + x: u a '\(5yU fl(v.)) + Xg'fa, 

+ X x t ^5{y^-fM) us :=J tus . (46) 

All the terms of this decomposition are (morally) of order three. Now, remember that 
we wish to tackle the case 37 > 1, so that it seems actually wise to invert S at this point, 
and we get 

f S tu dxi /,(„„) = Xf/fM + XZ>'{v, fM) s + X*:™(y*j ■ f>(y))s + A*.( J), (47) 

where J tus is given by (|46|) . Notice once again that we have obtained a decomposition 
valid for some smooth functions x and y, but this decomposition can now be extended 
to an irregular situation up to 7 > 1/3. 

In a natural way, we will use ( T471) as the definition of the integral in the prescribed 
context of a 7-Holder process with 7 > 1/3. To conclude this heuristic reasoning, let us 
summarize the different hypotheses we have (roughly) raised during the procedure: 

• The process x generates four operators X x , X ax , X xa et X xx , which satisfy the 
algebraic relations (13"6"I) . (T4"4l) and (HB]) . As for the Holder regularity of those 
operators, X x admits the same regularity as x, X xx twice the regularity of x, 
just as X ax and X xa (even if one must change the space one works with, according 
to the above point (a)). 

• The increments (Sy) ts can be decomposed as (5y) ts = X x s y x + yf s , where is 
twice more regular than y. Besides, according to (a) again, the process y must 
evolve in a space B ajP , with a > 0. These remarks will give birth to the spaces 
Q K 

• The functions <7j are regular enough (to be precised below). 
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Remark 4.3. If one has a look at the constructions established in [49], it seems more 
natural, at first sight, to search for a decomposition of the integral based on the (twisted) 
iterated integral 

Xts ,ij := f S tu dx\ B{X*J ® Id) = J S tu dx l u B (£ S uv dxi ® Id) , (48) 

rather than on the area Xf® we have introduced in ( 1431) . In a way, the definition of 
Xf x is actually more consistent with the general iteration scheme of the rough path 
procedure. Nevertheless, when it comes to applying the results to a ffim x (with Hurst 
index H G (1/3,1/2)) for instance, it seems difficult to justify the existence of the 
iterated integral (148]) . According to our computations, this difficulty is due to a lack 
of regularity for the term S uv in (1481) . Indeed, if one refers to [2], the definition of the 
integral would require a condition like 

N\S UV — S uu ; C(B a>p , 13 a , p)] ^6 W ~~ v \ i 

for some v > 0, but this kind of inequality cannot be satisfied in this general form, since 
the Holder property (114"]) of the semigroup requires a change of space. This is why we 
have turned to a formulation with Xff, which is made possible by the introduction of the 
operator X™ (defined by (1371) ) in the decomposition (l38~l) . As we shall see in Section El 
the definition and the estimation of the regularity of X xx are much simpler, since this 
can be done by means of an integration by parts argument. 

4.2. Definition of the integral. In this subsection, we will only make the previous 
assumptions and constructions more formal. From now on, we fix a coefficient 7 > 1/3, 
which (morally) represents the Holder regularity of the driving process x. The definition 
of the rough path above x associated to the heat equation is then the following: 

Hypothesis 3. We assume that the process x allows to define operators X x ' 1 , X ax ' 1 , 



X xa ' 1 , X xx ^ E {1, . . . , A^} ), such that, recalling our Notation \4. 2\ ' 
(HI) From an algebraic point of view: 

SX*'* = (49) 

X x,i = X ax,i + 5x i ( 5Q ) 

5X xa ' 1 = X xa '\a ® Id) + X x '\a <g> Id) (51) 

SX xx ' ij = X x '%5x j ). (52) 
(H2) From an analytical point of view: if lap > n, then 

X x * e (%(C(B P , B p )) n C2(£(B a , p , B a , p )) n Cr n/{2p \C{B p/2l B p )) (53) 

X ^eC] +a (C(B a ^B p )) (54) 

xxa ,i e C T+a-n/(2p) (£(B ^ x ^ ^ R (H^^ x B ^ ^ (55) 
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X xx,n e c?(£(B p , B p )) n C 2 ^(£(i3 Q , p , n C 2 27 (£(£ Q , P , B p )). (56) 
VFe mi/ denote by X = (X x , X ax , X xa , X xx ) the path so defined. X belongs to a 
product of operators spaces, denoted by CC" hK ' p , and furnished with a natural norm build 
with the norms of each space. 

The formal definition of controlled process takes the following form: 

Definition 4.4. For all a G (0, 1/2), k G (0, 1), we define 

Q*, P = Q K a ,M T]) = {ye #([0, T], B a , p ) : (5y)ts = X*fyf + yt 

y*'* G £*>([(), T],£ a , p ) r\C?([0,T\,B p ), yt G C 2 7 ([0, T], B a , p ) n C 2k ([0, T], B p )}. 
We will call Q^p the space of k- controlled processes of B a<p , together with the norm 

N 

A% QaJ = A%; C?(B a>p )} + {A'//-': Cl{B a>p )] + W; } 

i=l 

+ Ar[ytt;C 2 '(i3 a , p )]+Ar^;C 2 2K (i3 p )] ! 
where the time interval [0,T] omitted for sake of clarity. 

Observe that, in what follows, we will only consider the spaces Q^ p , with 2np > 1. 
We can now show how nonlinearities of the form given in Hypothesis [2] act on a 
controlled process. 

Lemma 4.5. Assume that /j G A" 2 /or i = 1,...,N and let k G (1/3,7). If V £ 2k P 
admits the decomposition 8y = X X)l y x,% + y^, then the increment 5fi(y) can be written as 

S(fi(y)h = K ® *0(v, //(y))- + • ■ /<(y))- + /i(v)L (57) 

with fi(y)^ = fi{y)^ 1 + fi(yY' 2 + fiiy)^ 3 , where the elements fi(y)^' k are given by [Sty) 
and ( [Top . Moreover, one has 

Af[f t (yf l ;C! K (B p/2 )] < c fiX {-A%; Ci(B a ,p)] 2 + M[y\ Q^] 2 } (58) 
Mmr, Cl"{B p )\ < c f , x Af[y; & j , A' ./',(// C 2k (B p )] < c /)X .A%; . (59) 

Proof. This refers to the decomposition fl38|) . The estimate of fi(y)^ 2 is obvious, while 
the estimate of fi(y)^ 3 stems from the hypothesis 054p . As for fi{yy ,x , notice that 

\\fi{y)uh p/2 < c f \\(5y) 2 ts \\t3 p/2 < \\(Sy) t s\\ 2 Bp < IK^tolllp + \\atsy s \\ 2 Bp , 
and the result then comes from the property (TT^j) . 

□ 

We are now in position to justify the use of (14 7p as a definition for the integral: 
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Proposition 4.6. Let y G Q«p([0, T]) admitting the decomposition 5y = X x,l y x ' 1 + y\ 
with k G (1/3,7) an d P G N* such that 7 — k > n/(2p). Assume that f = (fi, . . . , f^) 
with j\ G Xi for i = 1, . . . , N. We set, for all s < t, 

J ts (dx f(y)) = ,Y;';7,(//,) + X?% //(y)). + X^(y^ ■ //(„)). + A ts (J), (60) 
where we recall our Notation \4-l\ for f[, and with 

Jtus = XgtfMl) + X*:%{5y) us , //(„.)) + X^\y u , *(//(*))«) 

+ X^5(y^.fi(y)) us , (61) 



the term f{yy being defined by the decomposition (57|). Then one has: 

(1) J(dx f(y)) is well-defined and there exists z G Q£ )P ([0, T]) such that 5z is equal 
to the increment J{dx f(y))- Furthermore, for any < s < t < T, the integral 
Jts(dxf(y)) coincides with a Riemann type integral for two regular functions x 
and y. 

(2) The following estimation holds true 

M[ Z] Q:, p ([0,T])] < c f>x {l + N[yA{B K , p )f + T"M[y; Q^] 2 } , (62) 

for some a > 0. 

(3) For all s <t, 

j ts {dxf{ v ))= lim {^UMy^+ x tZt k (yJM)t k 

I— *u 

+ K x ?L(y x ' j >fl(y)k}, (63) 

where the limit is taken over partitions A[ S)t ] of the interval [s,t], as their mesh 
tends to 0. 

Proof. The fact that Jtsidx f{y)) coincides with a Riemann type integral for two regular 
functions x and y is just what has been derived at equation (HTI) . As far as the second 
claim of our proposition is concerned, it is a direct consequence of Hypotheses [2] and [3j 
together with the estimations of Lemma [4.51 Let us check for instance the regularity of 
J: 

• for X^ifiiyY), we get, by §3J and §HD, 

Af[x^u(y) 12 + Myf 3 y,cr 2K m < c/ , x .a% sy, 

while, owing to flo^l) and 



XiX*>JM l \c:t 2K - n/(2p \B p )} < c /iX {N[y-A^ P )f + M[y; frj] . 
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• for X xa ' l ((5y), //(?/)), the hypothesis (1331) gives 

• for X xa ' l (y,5(fl(y))), one has, by (1331) again, 

AA[X-(y,5(/;(y)));C3^ 2K ^ (2p) (i3 p )] 

< c /> x^[2/;C 1 (B K)P )K[ 2 /;C 1 K (^)] 

< c^A^C?^)] {AT[y;C 1 °(i3 K , p )] + A%; Qjy} . 

• for X xx ' ij S(y x ' j ■ fl(y)), we deduce from (135]) that 

< c /jX {AV^Cf^)] +Arr^C 1 °(i3 K , p )]Ar[y;C 1 K (i3 p )]} 

< c /lX 1 1 + A%; C°(i3 KjP )] 2 + A%; Q^] 2 } . 

Moreover, thanks to the algebraic relations stated in Hypothesis [3] and the decomposi- 
tion ( 1371) . it is easy to show that 

J = -5 {X x >U(y)) + X™>% fi(y)) + X x ^(y^ ■ //(y))) . 

Therefore, J G Ker<5 fl C^(B P ), with /i = 7 + 2k — n/(2p) > 1, and we are allowed to 
apply A. Besides, using the contraction property ( 122]) . we get 

Af[k(J);Cr 2K - n/i2p \B P )} < C/ ,x |l + A%;C°(£> KjP )] 2 +Af[y; Qlj) , 

and also 

Ar[A(J);Cr K - m/(2p) (i3 K , p )] < c />x {l +U[yA{B K , p )] 2 +U[y; Q^] 2 } . 

The regularity of the other terms of (1601) can be proved with similar arguments. As for 
the expression (163]) . it is a consequence of Proposition 12.11] since one can write 

J(dxf(y)) = (id -A£) (X'' i (f i (y)) + X^ i (yJi(y))+X m ^ ■ //(y))) . 

□ 

Once our integral for controlled processes is defined, the existence and uniqueness of 
a local solution for our equation is easily proved: 

Theorem 4.7. Assume that f = (f\, . . . , f^) with f± G X 3 for i = 1, . . . , N. For any 

pair (n,p) G (1/3, 7) x N such that 7 — n > n/(2p), there exists a time T > for which 
the system 

(5y) ts = J ts (dxf(y)) , y = ^eB p , (64) 
interpreted with Proposition \4 ■ 6] admits a unique solution y in Q* p ([0,T]). 
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Proof. This local solution is obtained via a standard fixed-point argument in the space 
of controlled processes. The procedure essentially leans on the estimation (1621) . The 
interested reader can refer to [28] for further details on the principle of the proof. 

□ 

5. Global solution under stronger regularity assumptions 

The aim of this section is to show that a regularization in the nonlinearity involved 
in our heat equation can yield a global solution. Specifically, this section is devoted to 
the proof of the existence and uniqueness of a global solution to the (slightly) modified 

system 

{h)ts= I Studx® S e fi(y u ) , yo = ip, (65) 



where f\ G X 3 , ip G B a , P for some a > to be precised, and e is a strictly positive fixed 
parameter. Owing to the regularizing effect of S e , we will see that such a system is much 
easier to handle than the original formulation (JMj). 

Note that we have chosen a regularization by S e in fl65|) . in order to be close to 
Teichmann's framework [36]. However, it will be clear from the considerations below 
that an extension to a convolutional nonlinearity of the form 

[ky)}(0= ! K&rifMv^dv, eel", 

with a smooth enough kernel K, is possible. The technical argument which enables to 
extend the local solution into a global one are taken from a previous work of two of the 
authors [TOl. 



5.1. Heuristic considerations. The regularizing property (fTBl) of the semigroup S e al- 
lows us to turn to a decomposition of J"* S tu dx^u S £ fi(y u ) similar to the finite-dimensional 
case, or otherwise stated written without the help of the mixed operator X xa . Indeed, 
let us go back to the decomposition fl38|) : 



S(fi(y))ts = (5x) ts y* ■ ft(y s ) + 



atsVs ■ f'iiVs) 



+ Vl ■ fiiVs) + {Xtryf) ■ fKVs) + / dr [f!(y s + r(6y) ts ) - f>{y s )} ■ (Sy) ts 



o 



(66) 



but this time, let us consider the whole term into brackets as a remainder term evolving 
in Bp (or maybe B p / 2 ), and denote it by fi(y)^. This point of view is for instance justified 
if we let the process y evolve in £> l p , insofar as, for any s, t G /, 

htsVs ■ fi{y S )\\B v ^ \t ~ S\ \\ oc\\y S \\ Bup < \t ~ S\ 2K lif" 2 * ||//||oo||!/«||Bi J> . 
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For obvious stability reasons, the strong assumption y s G Bi tP then implies that the 
residual term steming from the decomposition of J"* S tu dx$ S £ fi(y u ) should also be seen 
as an element of B\ )P . This is made possible through the action of S £ . Indeed, owing to 
( TT3|) . one has 

WSeifiyf)^ < ce" 1 \\f(yf\\ Bp , for some constant c > 0. 

5.2. Definition of the integral. According to the above considerations, only the pro- 
cesses X x ' 1 , X ax ' 1 and X xx ' 1 will come into play. Therefore, let us focus on the following 
simplified version of Hypothesis [3J 

Hypothesis 4. We assume that the process x allows to define operators X x,% , X ax,% , 



X xx m G {1, . . . , N} ), such that, recalling our Notation \4.S\ ' 
(HI) From an algebraic point of view: 

k 1 ' 8 = (67) 

X x,i = X ax,i + 6x i (gg) 

SX xx ' ij = X x >\5x j ). (69) 

The operators X x ' 1 and X xx,lj commute with S £ . (70) 
(H2) From an analytical point of view: 

X* G CJ(C(B P , Bp)) n C1(C(B 1>P , B 1>p )) n Cr n/(2p) (C(B p/2 , B p )) (71) 

X^ eC l 2 ^{L{B 1>p ,Bp)) (72) 

X x*M e c?(jC(B p , Bp)) n C 2 27 (£(B 1)P , B liP )). (73) 



Remark 5.1. The assumption (1701) is trivially met when x is a different iable process and 

X x ' 1 is defined by X^f = J* St u dxu ■ It will remain true in rough cases, following the 
constructions of Section O This commutativity property will be resorted to in the proofs 
of Propositions 15.31 and 15.41 

The notion of controlled processes which has been introduced in Definition 14.41 can 
also be simplified in this context: 

Definition 5.2. For any k < -y, let us define the space 

Q K , P ={ye CJ{B 1)P ) : (5y)ts = Xtftf* + yt, jf>* G C^{B 1)P ) n C«{B 1)P ), y* G C 2 2k (B 1>p )} , 
together with the seminorm 

With this notation, one has Af[y; Cj(Bi tP )} < c x J\f[y; Q K , P }- 

In the following two propositions, let us fix an interval I = [a, b] and denote |/| = b — a. 
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Proposition 5.3. Let y G Q KtP (I) with decomposition 5y = X x ' l y x ' % + y^ , for some 
(n,p) G (1/3,7) x ^* such that 7 — k > n/(2p) and initial value h = y a G &i, p - For any 
ip G Bi, p , define a process z by the two relations: z a = ip and for any s < t G I, 

(Sz) ts = j u (dx® s e fi( v .)) = xpSifM + x^sM' J ■ f!(v.)) 

+ A ts (X^SMyf + X''-»S, Si,/--' ■ f>(y))) , 

where fi(y)* stands for the term into brackets in UJh}) . Then: 

• z is well-defined as an element of Q K>P (I). 

• The following estimation holds: 

Af[z;Q K , p (I)\ < ce- 1 {l + |/| 2( ^ k) A%; Q K , P (/)] 2 + WHl^} , (74) 

for some constant c > 0. 

• For any s < t G I, (Sz) ts can also be written as 

(5z) ts = lim { X Z ltk SJi(yt k ) + *Z?LS e (y% j ■ fKyJ)} inB 1 , p . (75) 

\ V Mh°t k €V [e , t] 

Proof. Let us focus on the estimation of the residual term 

zt = X;:-"Si«:-< ■ f.(y.)) + A ts (X'-'S f,(yy + X-'-'S Ay'-' ■ f[{y))) . 
First, using (173"]) and (fI3|) . we get 

\\X^SM> 3 -ti(ys))W v < c x \t-s\^e- l \\y^-f>( ys )\\t3 P 

< c x \t-s\^e- l \\ y ^\\ B , p 

< c x \t-s\ 7rf e- 1 jsf\y;Q KiP {I)]. 

Secondly, write /,(y)« = f\{y)^ + /,(„)«.*, with f^s = a ts y s ■ ft(y s ) + y\ s ■ fl(y s ) + 
(X^Y/) ■ fi(y s ), fMf = Jo dr [fi(y s + r(Sy) ts ) - f{(y s )} ■ (6y) ts , and notice that 

Wx^sjMlilh^ 

< \t - e^wm^ 

< \t - e- 1 {\\(a us y s ) ■ //(y.)^ + \\{X«fyf) ■ fl(y s )\\ Bp + \\yt • //(y s )|| Bp } 

< |t - uFe- 1 {\u - s\ \\y s \\ Bhp + \u- s\ 1+ ^ \\yf\\ Bhp + \\yt\\ Bl J 

< \t ~ uf'e- 1 [\u - s\ 2K Af[y; Q K , p (I)} + \u- s\ {a%; Q KJp {I)\ + \\h\\ B ^)) 

< \t ~ sl^e" 1 {a%; Q R , P {I)\ + \\h\\ Bl , p } , 
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while, owing to 


(HOD, 


XtuSefiivfusllB! 


p 


< e-^t-ur 


n/(2p) 




n/(2p) 




n/(2p) 


< e- x \t-ur 


n/(2p) 




-n/(2p) 



\\fi(y)ush p/2 

IIO*2/)w||L 



i|2 I II ||2 
us\\B p \\ a u S ys\\l3 p 



\h\\k, P 



|J|2(l-7) || 









\u-s\ 2 ~<M[y;Q^ p (I)] 2 + \ U 
Even more simple estimations based on f!73|) give 

\\x^s e 8{y x>j ■ mUW P 

'l +A%Q K , P (/)] 2 + \lt K ^[y\Q.M\ ■ \\h\W P \ ■ 



< s' 1 It 



|27+k 



Thanks to the contraction property (|22|) . we now easily deduce 

M[z^Cl K {I)\ < ce- 1 {l + \I\^~ K) A%;Q K , P (/)] 2 + |/| 2(1_K) \\h\\ 2 BltP 

The estimation of Af[z x,i ; C? K (I; B hp )} can be established along the same lines. As for 
HJ), it is a consequence of (12. lip , together with the reformulation 



Sz = (Jd-k$)(X'* i S e f i (y) + X m >*S e (3r> ■ ti(y))). 



□ 



In order to settle an efficient fixed-point argument in this context, the following Lip- 
schitz relation is required: 

Proposition 5.4. If y,y G Q KjP {I) with y a = y a , and if we denote by z, z the two 

processes in Q K}P (I) such that 



z = z = y 



and Sz = J(dx® S e fi(y)) , Sz = J{dx® S £ f t (y)), 



then 



M[z - z- Q KiP {I)] < c x e~ l \ir K Af[y - y; Q K , P (I)] 



|2(7-/t) 



W[y;Q K M} 2 + *f[y,Q K M)} 2 } + |/| 2(1 " K) \\h\\j^ \ ■ (re) 
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Proof. One has, for any s,t E I, 

S(z - ~z) ts = x*fs e Ui{y s ) - fi(y s )) + X x t ^SM' j ■ fM - <C J ■ ti(y s )) 

+ A ts {X*>*SMy)l - /,(*/)«) + ■ f'M - y xJ ■ fl(y))) ■ 

Let us only focus on the more intricate term, that is to say X x,l S £ (fi{i))^' 2 — fi(y)^ 2 ) 
where, according to the notations of the proof of Proposition 15.31 

U{y)\s = I dr [fl(y s + r(Sy) ts ) - //(y.)] • (Sy) ts . 

Write 



o 



fM's - fiiyfts = [ dr [f'{y s + r(5y) ts ) - //(y.)] • 5(y - y) ts 
Jo 

+ (Sy)ts -S(y-y)ta- drr dr' f"(y s + rr'{8y) ts ) 
Jo Jo 

• 1 pi 



+ mi- drr dr'[f;{y s + rr'{8y) ts )-f;{y s + rr'{8y) ts )} 



In this way, 

\\fiiy)u - fi{y)tsh P/2 < \\S(y - y)tsh P {\\{Sy)tsh P + ll(<5£Uk} 

+ \\($y)ts\\ 2 Bp {\\y s - yshoc + \\yt - mh^} ■ 

Now 

\\8{y-y)tsh P ^ \\Hv - ifitshi* + \t~s\ \\(y s -y s ) - S sa (y a -y a )\\ Blp 
< \t-s\^N[y-mQ K , p {I)i 

while 



\\{5y) ts \\ Bp < \\(5y) ts \\ Blp + \\a ts (8y) sa \\ Bp + \\a ts S sa h\\ Bp 
< \t-s\ K \\ir K JV[y;Q K , p (I)] + \I\ 



l-K 



and finally 



hs-ysho. ^ \\y s - y s \W, p ^ \\y s -y s - s sa (y a -y a ) \\ Bli 

< \I\ 7 ~ K J\f[y — y; Q K , P (I)]. 
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This easily leads to 

*) M -/*(yf;^M < \ir K M[y-mQ K , P (i)} 

{l + |/| 2(7 " K) {A%; Q K , P (/)] 2 + A%Q K , P (/)] 2 } + |/| 2(1 " K) 11*11^} . 
Inequality ( l76i) now follows from standard computations based on Hypothesis HI 

□ 

We are now in position to state the expected global result: 

Theorem 5.5. Let fi G X3, for i G {1,...,N}. Under Hypothesis^ let (k,p) G 
(1/3, 7) x N* such that 7 — k > n/(2p). For any T > 0, for any ip G £>i, p , the differential 
system 

{6y)u = Ju{dx®S e fi(y)) , y = 4, 
interpreted with Proposition ^ '.3[ admits a unique global solution in Q KtP ([0,T]) . 

Proof. From the two estimations ( 1731) and (1761) . the patching argument is exactly the 
same as in [TOj, Theorem 4.16]. It consists in controlling both the norm of the initial 
value and the norm of the process as a controlled path on each successive intervals. For 
sake of conciseness, the reader is refered to the latter article for a detailed proof of the 
statement. □ 



6. Application 

As it was announced in the introduction, the goal here is to apply the previous abstract 
results of both Sections [3] and H] to a fractional non linearity given by the formula 

N 

X t (<p)(0 = Y,x l M^¥>(0), (77) 
i=i 

with a d- dimensional 7-Holder process x = (x 1 , . . . , x N ) with 7 > 1/3, and Oi some 
smooth elements of X2, as defined in Hypothesis EJ 

To this end, we know that it suffices to construct, from x, a path X = (X x , X ax , X xa , 
X xx ) which satisfies Hypothesis [31 Indeed, the latter assumption clearly covers Hypoth- 
esis [U of Section [31 

As usual in this paper, we shall proceed in two steps: we first work at a heuristic level, 
that is with smooth processes, and try to obtain an expression which can be extended 
to irregular situations. We then check directly Hypothesis [3] on the expression obtained 
in the heuristic step. 
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6.1. Heuristic considerations. Assume for the moment that x is a smooth K - valued 
function. Then the operators X x , X ax , X xa and X xx are defined by the formulae 

*2'(¥0(O= fs t M{t)d< , X?' i (<p)(t) = fatM&dxl (78) 

J s J s 

X t ?'W)(£) = [* S tu ((a us p) ■ m) dxi (79) 

J s 

x?r j (m) = f sum) k (so) 

Set now x\ s = f f dx u ® (&c) us . Then a straightforward integration by parts argument 
yields the following expression for the increments introduced above: 

= (5x% s + f AS tu (5x\ s du (81) 

J s 

X?r = j\s tu {W)usdu (82) 
X?r = [ Xg (AS*,® Id) du (83) 

J s 

= x 2 J j + f AS tu x^du. (84) 

J s 

These are the expressions that we are ready to extend to irregular processes. 
Let us only elaborate on how to get ( |83|) . Actually, it suffices to notice that 



t ft 
S t u{{ausV) ■ il>) dx l u = - / d u (X^)((a us <p) ■ ip), 

J s 

where, in the last integral, the partial derivative d u only applies to the operato X 
Then 

- J d u {X^){(a us (p) ■ ip) = [-X^ l {(a us (p) ■ tp)]] + / duX^(d u (a us cp) ■ ip) 

t 

xA 



x.i 
tu " 



duX^((AS us if) -ip). 

Remark 6.1. At this point, it is not clear that the integral expressions J"* AS tu (Sx l ) us du,... 
give rise to operators defined on <B QjP . For the moment, we only consider those expres- 
sions as operators acting on C^°. The extension to any space B a ^ p will stem from a 
continuity argument (see the proof of Proposition 16.21) . 
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6.2. Definition of the heat equation rough path. In a natural way, in order to 
extend expressions (1ST]) - (151)) to a Holder path x, one has to suppose that this process 
generates a standard rough path, that is to say: 

Hypothesis 5. We assume that x allows to construct a process x G C 2 27 (M n ®M n ) such 
that Sx 2 = 5x <S> 5x, or in other words 

(Sx 2 '% us = (6x% u (5xi) us , /../ 1 V. (85) 

This allows us to state the main result of the section: 

Proposition 6.2. Under Hypothesis^ the operators X x ' 1 , X ax ' 1 , X xa ' 1 , X xx ' lj defined 
by l[8l\) -(84\), can be extended to a path X which satisfies Hypothesis^ 

Proof. We have to check both the algebraic and analytic assumptions. 

Algebraic conditions. The verification of ( I49l) -(l52l) is a matter of elementary calculations. 
For instance, let us have a look at relation (1521) . For all s < u < t, one has 

(5X XX >% US = a£« _ aA« _ StuX ^ + f AS tv {x 2 ™ - x^) jdv. 



Then, by (J85J), this expression reduces to 



tus 

(Id -S tu )x 2 ^ + (6x% u (6xi) us + / AS tv (x 2 ^ + (6x% u (6xJ) us ) dv 

J u 

(5x l )tu+ / AS tv (d~x l ) vu dv (5x J ) us = X^(5x 3 ) us . 

J u 

Analytical conditions. Let us examine the regularity of each operator individually. 
Case of X x ' 1 . The norms at stake here are 

M[X x ' i \ C^(C(B P , B p ))} (86) 
N[X X <\C1{C{B^B K , P ))} (87) 

Af[X x ^Cr n/2p (C(B p/2 ,B p ))}. (88) 
In order to establish those regularity results, let us first rewrite (18T1) as 



X£ = S ts (8x% s - f AS tu (5x% u du. 

J s 



Then observe that (186j) and ( 1871) are obtained thanks to the same kind of arguments. 
We thus focus on (1871) for sake of conciseness. But the latter norm can be bounded easily 
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by noticing that: 

K%)lk, p < \\Sts&)\\B K J(5x% s \+ f \\AS t M\\ BK J(8x% u \du 

J s 

£ IMkJM| 7 (\t ~ S ! 7 + J \t ~ < IMkJMW* - s| 7 , 

which holds for all k > 0. Along the same lines, in order to prove ([58]) . we use the 
fact that \\S ts (ip)\\ Bp <\\ip\\B p/2 \t-s\- n / 2 P and that \\AS t3 (<p) \\ Bp < \\<p\\ Bp/2 \t - s]- 1 '^. 
Then we obtain 

ii^/Mik<ii^ik /2 ii^ii 7 |t- s r-^ 

for all p such that 7 — njlp > 0. Those estimations give the required bound fl88|) . 

Case 0/ X ax '*. We should now check that ( 15^1) is verified in our setting. To this aim, 
write X aX)l as 

x ts' 1 = (k s (5x l )ts - / AS tu (Sx l )tudu. 



Then 

ft 



IWfoOlk = KMIkl(<^%l + / US tu {y)\\ Bp \{8x*) tu \du 

J s 

and using the semigroup estimates 

iMvOik < iMkj*-*r iia^ u mii Bp < iMkjt-«r i+K 

we easily conclude that 

N[X ax >\ C(B K>P , B p )] < |* - s| 7+K , (89) 
which is the expected regularity result. 
Case of X xa ' 1 . Going back to (155]) . one must prove that the following norms are finite: 

tf [X xa,i. (%+«-»/<*>) (£(B Ka> x B p , Bp))], and A^r'; C 2 7 (£(B K , P x B K , P , B K , P ))]. (90) 
lo do so, write X ts as 

t/ s 

We deduce 

■A/Wfo V); Bp] < M[X^- cr n/(2p) (C(B p/2 , B p ))]M[{{a ts ^) • V); B p/2 ] 
+AT[X^ i ; C 2 7 - n/(2p) (£(B p/2 , Bp))] jT |u - s|W[((AW) • ^); B p/2 ]du 

where 

AT[((a 4 ^) ■ V); Bp/2] < A/K^l B P ]N{^ B p ] <\t- S | K A%; B K , P ]M[^] B p ] 



NON-LINEAR ROUGH HEAT EQUATIONS 35 

and 

M[((AS us <p) ■ j>); B p/2 ] < M[AS US ^ B P ]M[^ B p \ <\u- s\- 1+ "M[<p; B K , P }M[^; B P \. 
This allows to conclude that 

< Cr n/(2p) (C(B p/2 , B P ))W&; B K:P ]Af[ip; B p ] \t - 8 [r+™/to), 

and the first of the required bounds in fl90l) follows. For the second one, we have 

M[X£>\k V); B K J < N[X X ^ C1(£(B K , P , B K , p ))W[((a tsV >) ■ ^B K , P ] 

+ N[X X *-C!(C(B K , P) B K , P ))} ! \u- S yAf[({AS us <p)-i>);B KtP ]du, 

J s 

and using the algebra property of B K>P , we get 

Af[((a ts tp) ■ ^);B KjP ] < A%; B KiP ]Af[i>- B KjP ] 

and 

N[((AS us <p) ■ rj>); B K J < N[AS us <p; B^ P ]M[^ B KtP ] <\u- s^Af^; B K , P ]N[^ B KtP ] 
so that 

M[X^\^)-B P ] 

< M[X X '\ CJ(C(B K:P , B K , p ))]M[ip; B K , p ]M[t, B KtP ]{\t - s\t + [ \u- s^du). 

J s 

The second estimate follows. 

Case of X xx,l i. We must estimate the norm 

U[X^;Cl\C{B p ,B p ))i (91) 

and also Af[X xx ' ij ; C%> (£(B a , p , B^))} and J\f[X xx > ij ; C% \C{B aiP , B p ))}. We focus on (|9"T]). 
the others terms having similar behavior using the algebra property of B 0)P and the 
Sobolev embedding B a , P C B^. 
First, write X xx,v as 

x xx,v = s tsX 2,v _ J Ag tu \r*w + (5x l ) tu (5x 3 ) us ] du. 
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From this expression, we immediately get 
M[X^^)-B P } 

< c x ^f[S ts (ip); B p ]\\t - s\ 2 ^ + J MIASM; B p ][\t - u\ 2 ^ + \t - u|> - s\"*\]du 

< c x ^M[tp;B p ]\\t- s| 27 + N[p\B p ] J \t-u\- x [\t-u\^ + \t-u\i\u-8\i\]du^ 

< c x M[<p;B p }\t-s\ 2 \ 

This gives the expected conclusion J\f[X xx > ij ; C^(C(B P , B p ))] < oo. 



□ 



We are thus in position to apply the abstract Theorems 13.61 and 14.71 in order to solve 
the heat equation for a general rough path above x: 

Theorem 6.3. Let x = (x^\ . . . ,x^) a d- dimensional ^-Holder path fj > 1/3) satis- 
fying the rough path hypothesis and consider the infinite- dimensional noise X build 
on x through the formula (f77[ ). If f G C 3,6 (1R; M), then the stochastic differential system 

{Sy)ts = J S tu dX u (y u ), y = tp e B KjP , (92) 

interpreted with Proposition 1 3. S\ if 7 > 1/2 and Proposition \4 . 6] if 7 G (1/3, 1/2], admits: 

• A unique global solution in C°' K ([0, T], B KtP ) if H > 1/2, where the pair (n,p) G 
(0, 7) x N* is such that H + k > 1 and 2np > 1. 

• A unique local solution in Q£ p ([0,T*]) if 7 G (1/3,1/2], where T* is a strictly 
positive random time and the pair (n,p) G (1/3,7) x is such that H — n > 
l/(2p). 

Remark 6.4. It is a well-known fact that one can construct a rough path (in the sense 
of Hypothesis [5]) above a ^-dimensional fractional Brownian motion B with Hurst pa- 
rameter H > 1/3 (see e.g. [SJ Q21 ESI SO])- This means that we can solve the heat 
equation (19"2"1) driven by this kind of process. 



7. Rough case of order 3 

To conclude with, and also to reinforce the feeling that our approach to the problem (jSj) 
is viable, let us say a few words about the case of a 7-Holder noise x, with 7 G (1/4, 1/3] 
only. We will not present the construction of the integral with as many details as in the 
previous section, and will stick to the broad lines of the calculations. 
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7.1. Construction of the integral. Fix an index 7 G (1/4, 1/3] which represents the 
regularity of x. In order to be allowed to invert 5 via Theorem 12.101 one must look 
here for a term of order 4, or more exactly of order 7 + 3k., where k is such that k < 7 
and 7 + 3k; > 1. The crucial point of the following construction lies in the (obvious) 
existence of a coefficient k G (0, 1/4) such that 7 + 3k > 1. Since k < 1/4, we can resort 
to the space £> 2ftiP and envisage the possibility of a solution evolving in this space. In this 
context, the operator X xa ' 1 that we have introduced in the previous section, and which 
was formally defined as X™' 1 = J* S tu dx l u (a us ® Id), becomes an order-three operator: 

X xa ^eC^ 2K (B 2 ^ p xB p ,B p ). 

Taking this observation into account, it seems then quite appropriate to consider the 
following space of controlled processes: 

Qk P = {y ■ $W ■ (k)u = x?/yf + Xu' ij v7> ij + vL 

{5y*)u = (Mh ■ yT' ji + Vu'\ 

y x > 1 e C^B^) n CZ(B P ) , y xx > ij G C^B^) n Cf(B p ) , y« G C 2 7 (i3 2ft , p ) n , 

together with its natural norm. One should notice the additional relation between y x 
and y xx which appears in the definition above, with respect to the rough case of order 2. 
This is reminiscent of the nilpotent algebra structure of [13], and also of the algebraic 
structures introduced in [TS1 [13, EZ] . 

According to our usual way to construct rough integrals, for the time being, x is 
assumed to be differentiate and the operators X x ' 1 and X xx ^ are defined by the formulae 

X t g = / Stu dx u , X ts 3 = / St u dx^ u [5x ) us . 

J s J s 

Besides, as in the previous section, the integer p is picked such that 4np > n and in this 
way, £>2 K , P becomes a Banach algebra. We shall then expand our integrals so that they 
can be extended to the case of an irregular noise. 

Assume that fi G A3, where A3 is defined at Hypothesis [2], and similarly to Nota- 
tion 14. 1[ set [//'(</?)](£) = V|cTj(£, </?(£)). If y G Q^p, an elementary Taylor expansion of 
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order 3 yields 

(Sf^U = (X:>iy? j ) ■ fM + (X^ k y^ k ) • #(y.) + a us y s • £(y s ) + \{8yf us ■ /f (y s ) 



1 rl 

2 

us 



+ Vis ■ fiiVs) + / drr dr' \f»(y. + rr'{8y) us ) - (y.)] ■ (6y) 
Jo Jo 



(5x% s ■ & ■ f'(y s ) + x 2 J^ ■ y^' jk ■ f((y s ) + a us y s ■ f((y s ) + l -{8yf us ■ f»(y.) 
+ ■ fi(v.) + X^ k {yT' jh ) ■ f&y.) 



1 r l 

2 

us 



+ Vis ■ fiiVs) + / drr dr' \f»(y. + rr\8y) us ) - £'(y s )] ■ (6y) 
Jo Jo 



and thus 



(Sfi(y))us = (Sx% s • //;•' • f'(y s ) + x 2 ^ ■ y^ k ■ £(y s ) + a us y s ■ f!(y s ) 

+ \{M) U8 yf j ■ (» r k )„,,r: J - ■ J7(v,) + mi. 

= (8xi) us ■ & ■ f[{y.) + x 2 J^ • {//;:•'•''" • //(y.) + y^ ■ y^ k ■ f>'(y s )} + a us y s • //(y.) 

+ ( 93 ) 

where fi{y)\ s is a residual term which is long (but easy) to write explicitly, and which 
can be estimated as 

Wt(y)S; *y < I* - *I 3K (i + A% Q 2 K J 3 ) • (94) 

In the calculation that leads to (1931 . we have introduced the operators 
x ts' 1 = / o tu c&4 , x^ x,1J = / a tu dx 3 u {8x % ) us , 

t/S J s 

while the notation a; 2 /- 7 stands for the usual Levy area x 2 ^ = f t dx J u (5x l ) us . As for the 
estimation (194jl . it is obtained by means of the following natural hypotheses: 

X a X ,i £ C 7+2« (£(B2kj} x ^ Bp) ^ and ^-.ii e C 2 27+2K (£(B 2lt , P> B p )). (95) 

Now, in order to be able to define Jt s {dx /(y)) for irregular processes, inject expression 
into the decomposition 

S tu dx u f(y u ) = X%\fi(y 8 )) + / S tu dx l u (8fi(y)) us , 
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to deduce 

' S tu dx u f(y u ) = .V;;'(. /•,(//,)) + X:r j (yf ■ f' 3 (y s )) + X?\y„ f t {y s )) 

+ ^r ijh (<rr j ■ I'M + yf ■ y x s ] ■ fttv.)) + r*, (96) 

with r ts = J S tu dx l u (fi(y)u S )- The operator X xxx ^ k is here defined by 

J s 

and we associate to this operator the (reasonable) regularity assumption 

Xxxxvk e C ?(£(B 2k , p , B 2k , p )) n C?(£(B P , Bp)). (97) 



According to the considerations of Section 14.11 it only remains to establish that 5r is 
a term of order 4, as a process with values in B p . Actually, we are going to show that 
5r G Cj +3h " n ^ p (B p ). It will then suffice to pick p large enough, so that 7 + 3k — n/p > 1. 

To compute 5r, one must assume the set of algebraic hypotheses (HI) in Hypothesis 
[31 and a further algebraic relation for X xxx : 

(<)X 11 ■■'■>'■),,<, = + X^ k (5x% s . (98) 

As far as the regularity assumptions are concerned, on top of condition ( 197]) . we have to 
modify a little the set (H2) in Hypothesis [31 which becomes: 

X x > 1 G C2(£(B P , B p )) n C 2 7 (£(B 2k , p , B 2KiP )) (99) 

j^, G Cr n/(2 p)(£(gp/2) Bp)) n C 7-n/ P(£(Sp/3) ^ )} (10Q) 

g C 2 27 (£(B P , B p )) n C 27 (£(B 2ft , p , B 2(tjP )) n C 27 - n/(2p) (£(B p/2 , B p )). (101) 

Besides, in order to clarify the presentation of this developpement, we will have recourse 
to the notation pa to signify "congruent to a term of order at least 7 + 3k — n/p" , or 
in other words: for all h, I G C3, h ~ / <^=> /i — / G Cg +3K n ^ p (B p ). Recall that, with this 
convention, our aim is to establish that 5r pa 0. 
Going back to (l9"6"j) . one has, by (1201) . 

-(5r) tus 

= • S(f(y))us) + (5X**'% us (yf • #„.)) - A"/,;' '' / fV(//' ' . /;(„))„. 

+0X* ae ^*) tu .(^ • tf(y.) + tf 4 • • (j,.)) 
-XSr^Stf** ■ f k {y) + ^ • y x > j ■ f'M) us 
+(5X™>% s (ys, fM) - X^((Sy) us , f'{y s )) - X^(y u , S(fl(y)) us ). 
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One can already notice that the fifth and seventh terms of the last sum have the expected 
regularity. Thanks to the above algebraic relations, together with the decomposition 
( 193|) . we then deduce 

-(Sr) tus 

« • f'M)) ~ •^/'•'• , ( </;:•' • *i ■ Kf{v))us) 

+x^ k ((Sx% s ■ {yT' ij ■ fM + yf ■ <C J ■ /;'(//-)}) 

+X™'\a us y s , fl(y s )) - X^'\(Sy) us , ft(y s )) 

« v'r (<>■::■' ■ °j ■ fXvs)) - x^aku m)) 

-X?u ij {y x u l ■ 5(f'M)us - (6x k ) us ■ y x / ■ yf ■ //(y.)) 
« (yf ■ WjMU - (Sx% s ■ y x ' k ■ f"(y s )}) . 

It is finally easy to see that (Sfl(y)) us = (Sx k ) us y x ' k ■ f['(y s ) + fl(y)i s , where f[(yf is a 
term such that 

WML Bp/2] <\t-s\ 2K B[ y] Q K 2 J 2 , 

By means of Hypothesis (jlOip . this statement enables to conclude <5r pa 0. 

Let us turn now to our main aim, which is an extension of the integral to Holder 
processes with Holder exponents greater than 1/4. We first formalize the assumption 
on X into: 

Hypothesis 6. We assume that the process x gives birth to operators X x,t , X ax ' 1 , X xx > l i , 
X^'', X axx ^ , X xxx,t i k , for which the algebraic conditions (HI) in Hypothesis^ and the 
analytical conditions {U ty) - fT0l\) . I[95\) . (9l\ ), are satisfied, for some triplet (7, K,p) such 
that 

7 G (1/4, 1/3] , k G (0, 1/4) , p G N* , Ak P > n , 7 + 3k - n/p > 1. 

Just as in Section HI this hypothesis allows to give a sense to the rough integral at 
stake here: 

Proposition 7.1. Under Hypothesis^ and assuming that f — (f\, . . . , f^) with f\ G A3 
for i = 1, . . . , N, we set, for any y G Q« )P ([0, T]), J(dx f(y)) = (Id— AS) (J), where, for 
all s < t, 

Its = XZVM) + X?*> ij (yf ■ fjiVs)) + X?r(y s , f[(y s )) 

+ x*r ijk (y7> ij ■ f'M + //.:•' ■ • fk(y s )). (102) 

Then one has: 

(1) J(dx f(y)) is well-defined and there exists z G Q^pQO'^l) such that 5z is equal 
to the increment J(dx f(y)). Furthermore, for any < s < t < T, the integral 
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Jts(dx f(y)) coincides with a Riemann type integral for two regular functions x 
and y. 

(2) The following estimation holds true 

A/>; Q 2 \ P ([0,T])] < CftX {l + U[yA{^ P )f + T a M[y; Q 2 K J 3 } , (103) 

for some a > 0. 

(3) For all s < t, J ts {dxf{y)) = ]im\ A[s ^ Y,(t k )eA [s>t] J t k +it k - 



7.2. Resolution of the differential system. As for Theorem 14.71 the (local) resolu- 
tion of the noisy heat equation of roughness order 3 stems from a standard fixed-point 
argument based on the estimation fl 1 U3j) : 

Theorem 7.2. Under Hypothesis® and assuming that f = (f\, . . . , f^) with fi e 
i = 1, . . . , N, there exists a times T > for which the system 

(5y) ts = Jts(dxf(y)) , y = e B p , (104) 

interpreted with Proposition \ 7. 1\ admits a unique solution y in Q,2 Kp {[0, T]). 



This abstract theorem can then be applied as in Section El thanks to the construction 
of a path 

j^ax,i j^xx,ij j£xa,i j^axx,ij j^xxx,ijk^ 

from a 7-Holder process x = (x (1) , . . . , x {d) ) , with 7 > 1/4. The rough path Hypothesis 
[5] must simply be enhanced in: 



Hypothesis 7. Assume that the path x allows to construct two processes x e C 2 
W 1 ), x 3 G C 2 ^(R n ® K n ® W 1 ) such that 

5x 2 = 5x ® 5x , x 2 + (x 2 )* = Sx Cg> 



27/ 



5x 3 = a; 2 <S> 5x + 5x (8> £ 2 . 

The path X can then be rigourously defined via the transformations fl8Tj) -f l84"j) . together 
with the additional expressions 

X?r ij = f AS tu x^ du, 



x xxx, ljk = x s,m + y * AStuX 3,m du _ 
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